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ABSTRACT 

Recently several type lb supernovae (SNe; with the prototypical SN 2005E) have been shown to 
have atypical properties. These SNe are faint (absolute peak magnitude of ~ —15) and fast SNe 
that show unique composition. They are inferred to have low ejecta mass (a few tenths of a solar 
mass) and to be highly enriched in calcium, but poor in silicon elements and nickel. These SNe were 
therefore suggested to belong to a new class of calcium-rich faint SNe explosions. Their properties 
were proposed to be the result of helium detonations that may occur on helium accreting white dwarfs. 
In this paper we theoretically study the scenario of helium detonations, and focus on the results of 
detonations in accreted helium layers on low mass carbon-oxygen (CO) cores. We present new results 
from one dimensional simulations of such explosions, including their light curves and spectra. We 
find that when the density of the helium layer is low enough the helium detonation produces large 
amounts of intermediate elements, such as calcium and titanium, together with a large amount of 
unburnt helium. Our results suggest that the properties of calcium-rich faint SNe could indeed be 
consistent with the helium-detonation scenario on small CO cores. Above a certain density (larger 
CO cores) the detonation leaves mainly 56 Ni and unburnt helium, and the predicted spectrum will 
unlikely fit the unique features of this class of SNe. Finally, none of our studied models reproduces 
the bright, fast evolving light curves of another type of peculiar SNe suggested to originate in helium 
detonations (SNe 1885 A, 1939B and 2002bj). 
Subject headings: nucleosynthesis, hydrodynamics, supernovae 



1. INTRODUCTION 

Recently, a new type of peculiar type lb super- 
nova (SN) has been discovered (with the prototypical 
SN 20 05E and a full sample of 8 SNe; iPerets et all 
l2010al) . These helium rich SNe show several intriguing 
features. They show fast and faint light curves (B-band 
peak luminosity of -15) and dominant lines of calcium 
in the nebular spectrum. Only a small fraction of ra- 
dioactive nickel is found in their ejecta, and Si group 
elements seem to be completely absent in their nebu- 
lar spectra. In addition, the enviro nment of these SNe 
is typically old (jPerets et al.ll20ldaH bT); inconsistent with 
massive young progenitors, usually thought to be asso- 
ciated with type lb SNe. It was therefore proposed that 
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these objects may result from helium detonations that 
occur on helium accreting white dwarfs. Energetically, 
the moderate observed velocities of 11,000 km see -1 to- 
gether with the estimated ejected mass (a few tenths of 
a solar mass), are consistent with the binding energy of 
0.2 M Q of helium. 

Other peculiar low mass and even faster evolv- 
ing SNe, but bright (peak luminosity of ~ —18.5) 
have also been studi e d (Ide Vaucouleurs fc Corwinl 
1981 IChevalier fc Plaitl 119881 IPoznanski et all I20T0F 



Perets et al J2010d . These SNe have also been suggested 
to result from helium detonations. Although we shall re- 
mark on these SNe, their light curves are brighter and 
faster than those produced in our simulations, and are 
unlikely to be produced by these scenarios. 

In the past decades, several authors discussed the 
scenario of detonations in helium layers, accreted on 
carbon-oxygen cores, suggesting "peculiar" SNe with 
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curves (see e.g. IWooslev et al. 1 


1980; 


Nomoto 19801 119821 iWooslev et al. 1119861 Livne 


1990 


Livne & Glasner 


199ft 119911 IWooslev & Weaver 


1994 


Livne & Arnett 1 


19951). There are several main 


issues 



here which deserve close investigation. The first issue 
is the explosion mechanism and especially the question 
whether a detonation which develops first in the helium 
layer could ignite a second (successive) detonation in the 
CO core. The second important issue includes the nu- 
cleosynthesis and its impact on the shape of both the 
light curve and spectrum. The third issue is the evolu- 
tion of helium accreting white dwarfs to thermonuclear 
runaway. 

The question whether a second detonation occurs has 
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been discussed many times. If this happens, the process 
will lead to the disruption of the entire star, and the 
consequences will be completely different from those ob- 
tained in th e case where the core is not burning. Livne 
& Glasner (|Livne 1 [19901 ILivne fc GlasneFl [19901 fl~99l 
have studied this, where a small reaction network of 13 
species has been used. Under the assumption of spheri- 
cal symmetry, the CO core always experiences ignition at 
the center, after a converging pressure wave, which prop- 
agat es from the core-helium boundary, reac hes the cen- 
ter (jLivne 1 [l990t IWooslev fc Weaver1ll99l . However, 
in a more realistic scenario, helium detonation is un- 
likely to start at a spherical shell simultaneously. Rather 
it is more likely to ignite at a small area near the in- 
terface, which can be approximated by a point ignition. 
The problem becomes a two dimensional problem with 
cylindrical symmetry around an axis defined by the line 
connecting the ce n ter of the s tar and the ignition point. 
ILivne fc Glasner I (| 19901 . 11991D have shown two ways by 
which the core can be ignited in this case. Depending 
on actual parameters (mainly the density of the fuels 
near the interface), the sliding helium detonation emits 
a strong oblique shock into the adjacent CO core, which 
in some cases is strong enough to drive a CO detonation 
near the interface. Otherwise, converging shock waves 
propagate inwards and eventually converge on the axis, 
but in this case off center. Consequently, driving a sec- 
ond CO detonation from that convergence point is also 
probable if the density there is not too low. Note how- 
ever that those results were obtained for rather massive 
CO cores of 0.85 Mq and above. 

The main study of the second important issue, namely 
the role of nucleosynthesis i n the explosion and its ap - 
pearance, has been done by iWooslev fc Weaver I (|1994f ). 
who performed a series of one dimensional simulations 
using a detailed reaction network. They used CO cores 
between 0.6 Mq and 0.9 Mq with helium mantles be- 
tween 0.13 Mq and 0.63 Mq (where the total was k ept be- 
low th e Chandrasekhar limit). In agreement with ILivne I 
(|1990f ). they find that the core is being ignited in all 
cases by converging pressure waves. A common impor- 
tant observable feature of those models is their very fast 
light curve, which rises to maximum over roughly 10-12 
days. The amounts of 56 Ni vary among models, between 
0.2 Mq and 0.98 Mq, with strong correlation between the 
core mass and the amount of 56 Ni produced. Moreover, 
significant amounts of calcium and titanium isotopes are 
being synthesized. Interes tingly, similar results w ere ob- 
tained by 2D simulations (jLivne fc Arnett 1119951 ). where 
similar models were simulated under the assumption of 
point ignition. A more comprehensiv e study in both 2D 
and 3D was recently carried out by iFink et al.l (|2007l ). 
Their results confirm the main conclusions of the previ- 
ous multi-dimensional studies. 

Finally, the evolution of helium accreting white dwarfs 
to thermonuclear runaway is a very complicated subject. 
Moreover, the onset of helium detonation at the base 
of such helium layers is a speculative process, yet to be 
explored. Some progres s has been recently reported by 
IShen fc Bildsten I (|2009[) . In the context of our study, 
their main result is the minimal helium mass required 
for dynamical burning as function of core mass (Fig. 5 
there). They also point out that the exact composition 
of the accreted helium plays an important role in the 



runaw ay. Previous results, published bv lYoon fc Langer I 
<|2004f ) are consistent with the above. However, they fo- 
cused on the possible effects of rotation and found that 
rotation may inhibit the runaway. 

In this paper we ignore most of the above complica - 
tions and repeat the work of Wooslc v fc Weaver! (|1994[ ). 
but extend the range of core masses to lower values. Con- 
trary to previous work, which focused on the higher edge 
of the mass range, we are interested here in low mass 
progenitors which presumably can reproduce the obser- 
vations of SN 2005E-like objects. Note that part of the 
phase space of low mass cores was al s o rece ntly (and in- 
dependently) studied bv IShen et al.l (|201Q[ ). as we shall 
briefly discuss later. 

2. EXPLOSION MODELS 

2.1. Tools and initial configurations 

We use the hydro code Vulcan/ID (V1D) which con- 
sists of an explicit Lagrangian scheme. The reaction net- 
work is based upon REACLIB with 160 elements be- 
tween hydrogen and Ni, including neutrons, 1,2 H, 3 ' 4 Hc, 
7 Li, 7 Be, 8 B, 12 " 14 C, 13 " 15 N, 14 " 18 0, 17 " 19 F, 18 " 22 Ne, 
20 " 23 Na, 21 " 26 Mg, 23 " 27 A1, 24 " 30 Si, 27 " 31 P, 29 " 34 S, 

31-37Q 33-38^ r 36-39j£ 37_44 Ca 40 ~ 45 Sc 41_50 Ti 
44-51y ' 45-54Q. ' 48-55jyj n 49-58p e ' 50-59^ 53-64j^- 

(jWeiss fc Truran I 11990(1 . The initial models are hy- 
drostatic configurations, consisting of CO cores, with 
masses of 0.45, 0.5 and 0.6 Mq, and He layers of 0.15, 
0.2 and 0.3 Mq. The CO core consists of pure, equal 
mass fractions of 12 C and 16 0, while the He layer is 
pure 4 He (except for the model discussed in 12.31 which 
had 12 C mixed into the He layer). The CO core is 
assumed to be isentropic, with central temperature of 
10 7 K and density that fits the given core mass. The he- 
lium laye r is also isentropic, with bottom temperature o f 
2xl0 8 K (|Yoon fc Langer ll20MIShen fc Bildsten II2009D . 
Our temperatures a r e a b it higher than those used in 
IWooslev fc Weaverl (|1994f ) and therefore our densities 
arc slightly lower. The models are well resolved with 
more than 1000 spatial zones, where the CO core is repre- 
sented by sa 500 zones, gradually growing from 0.001 Mq 
at the center to 0.005 Mq at 0.05 Mq below the edge of 
the core, then gradually decreasing to 10 -4 Mq and re- 
maining at this value throughout the He layer. In TableQ] 
we detail the hydrostatic parameters of the cases we have 
simulated. 

The detonation at the base of the helium layer is driven 
artificially by giving 20 zones large positive radial veloc- 
ity of 10 9 cms" 1 . Usually this drives a detonation imme- 
diately, but, at very low density the detonation may die 
out after a short period of time. This suggests that the 
question of how/when helium detonation may be ignited 
spontaneously, should be studied separately under appro- 
priate conditions and smaller scales. In any case, when 
formed, the self sustained detonation is rather weak, 
leading to incomplete helium burning over most of the 
layer. As mentioned earlier, spherically symmetric simu- 
lations lead always to CO ignition at the center. To avoid 
this complication we artificially inhibit here the burning 
in the core. This is temporarily justified by the fact that 
in the 2D case the convergence of those waves is off cen- 
ter, and may not ignite the core when the density at the 
convergence region is low enough [see ILivne fc Glasner"! 
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(fl990l fl99lh for details]. 

The simulation is run without rezoning until the shock 
wave approaches the center (typically at ~ 0.1 Mq from 
the center, ~ 1 s after the detonation is ignited, sec 
Fig. [T]). From then on, as the He layer draws away from 
the core, the core is rezoned more and more coarsely. 
This is done in order to prevent the shocks, that are 
going back and forth through the core, from unnecessar- 
ily diminishing the time step (this does not affect the 
dynamics of the He layer). The simulation is run until 
10 5 s, at which epoch the densities in the He layer are 
low enough so that radiative transfer calculations can be 
carried out. Nuclear reactions are followed only above a 
temperature of 10 7 K, whereas weak interaction decays 
are followed throughout the simulation for all zones in 
the helium layer. 

2.2. Results 

The evolution of the density, temperature, and velocity 
profiles of a typical model C0.45Hc.2, having a CO core 
of 0.45 Mq and He layer of 0.2 M , is shown in Fig. Q] 
The initial velocity of 10 9 cms -1 at t = injected at the 
base of the He layer (in order to initialize the detonation) 
can be seen. At t = 0.1 s, a detonation front has already 
formed, the maximum temperature in the He layer is 
1.8xl0 9 K. Ati= 1 s, the outward going shock has swept 
through all of the WD, while the inward going shock is 
approaching the center. The maximum temperature in 
the He layer at this stage is 1.5 x 10 9 K, whereas in the 
core it reaches 2.5 x 10 s K. It is worthwhile to mention 
that when the shock reaches the center, the temperature 
there sharply rises to ~ 1.3 x 10 9 K. At t = 10 s, the ve- 
locity in the He layer is nearly homologous, and nuclear 
burning has almost ceased, as the maximum tempera- 
ture in the He layer has dropped to 1.5 x 10 8 K. From 
then on, as can be seen for t = 100 s and t = 10 5 s, the 
He layer continues expanding homologously, nuclear re- 
actions being exclusively weak-interaction decays. 

The compositions of our models are plotted in Fig. [2l 
The plotted epoch is 10 5 s, at which point the model is 
transferred to the radiative transfer calculation. For each 
model the figure shows (left panel) the composition by 
elements, and (right panel) the summed up composition 
of the species constituting the decay chains considered in 
the radiative transfer calculations (see Sec. [3]). 

The main nuclcosynthetic products from our simula- 
tions are summarized in Table [5] The main result of 
our simulations is the sharp drop of radioactive 
nickel products for cores below 0.6 Mq. Accord- 
ingly, larger amounts of 40 Ca, 44 Ti and 48 Cr are pro- 
duced, with almost no silicon group elements. 

In addition we note, that the remnant CO core is ex- 
pected to undergo pulsations due to the inward moving 
shock wave initiated by the detonation being repeatedly 
reflected from the center and the outer boundary of the 
core, but this phenomenon is not adequately resolved in 
our present simulations. 

2.3. Effect of initial composition on nucleosynthesis 

A major drawback of pure helium detonations is the 
fact that usually the burning goes all the way to 56 Ni, 
and there is no evidence for intermediate mass elements 
in the emerging spectrum. One way to overcome this 



drawback, as we demonstrate here, is to ignite the deto- 
nation at low densities. Once a pure helium detonation is 
ignited at densities below approximately 5 x 10 5 gcm -3 , 
the burning rate is low enough so that the composition, 
for the appropriate dynamic timescale, includes many 
intermediate mass elements and almost no iron group 
elements (Table Fig.©. 

Another possible way to generate a composition of 
ejecta dominated by intermediate mass elements is to 
ignite the detonation in a mixture of helium and carbon 
(or CO). For the relevant temperatures we obtain, the 
triple a reaction is slow compared to a capture on car- 
bon. Therefore, any pre existing 12 C will tend to capture 
the free a particles. It follows that, if there are sufficient 
carbon nuclei present at the onset of burning, the end 
products can be calculated by a straight forward argu- 
ment. 

The number N of available a particles per 12 C is given 
by: 

N = (X He /4)/(X c /l2) = 3(X He /X c ) (1) 

where Xn e and Xq are the mass fractions of He and 
12 C respectively. 

Initially Xc + Xjj e = 1: therefore the number N of 
available a particles a 12 C nucleus will ultimately accu- 
mulate is N = 3(1 - X c )/X c . 

We assume here that a captures simply continue up 
the a capture chain from 12 C. 

In this way, for any initial pre existing carbon abun- 
dance Ac, we can compute N and predict the end prod- 
uct. For initial carbon abundance Xc, the predicted 
major final nucleus in the aftermath of the detonation 
is given in Table [3] Once the final nucleus is produced, 
there are no more free a particles available for capture on 
heavy nuclei or for consumption by the triple a reaction. 

For a model in which we artificially assumed an initial 
carbon abundance of Xc = 0.3, there are seven free a 
particles for each carbon atom and indeed the end prod- 
uct of the detonation is found to be 40 Ca, in accordance 
with the predictions given in Table |3l 

A crucial issue for such a scenario is the origin of the 
pre existing carbon. Assuming that the outburst occurs 
on a carbon WD that accretes helium, carbon enrichment 
can take place if there is dredge up mixing at the bottom 
of the envelope prior to the ignition of the detonation. 
Preliminary ID models show that the helium envelope 
is indeed unstable to convection about a day before the 
runaway. We intend to examine this interesting possi- 
bility for the mixing process by performing ID and 2D 
simulations of the pre runaway evolution. 

3. RADIATIVE TRANSFER MODELS 

To connect the theoretical predictions from the ex- 
plosion models to the observed properties of SN 2005E, 
radiative transfer models which describe the generation 
and propagation of light within the ejecta are needed. 
Here we investigate different aspects. First we discuss 
models of the bolometric light curve describing the tem- 
poral evolution of the total light emission in ultravio- 
let, optical and infrared (uvoir) wavelength bands. Sec- 
ondly, we computed a sequence of synthetic spectra for 
the early epochs from the time of explosion to about a 
week after maximum light. 
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Fig. 1. — Evolution of density (solid red line), temperature (long-dashed green line) and velocity (short-dashed blue line) profile of model 
C0.45He.2, having a CO core of 0.45 Mq and He layer of 0.2 Mq. From t = 10 s, the velocity is homologous, as shown against the radius 
(dotted purple line). 
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3.1. Bolometric light curve models 

The code we use to predict UVOIR light curves from the 
explosion models provides a Monte-Carlo-based, time- 
dependent description of the propagation of radiative en- 
ergy packages through the ejecta. The procedure consists 
of a simulation of the transport of energetic 7-photons re- 
leased in the radioactive decays within the ejecta and the 
transport of optical photons generated after the interac- 
tion of the 7-photons with the matter. The procedure 
foll ows the description g iven in ICappellaro et aTl (|1997f ) 
and lMazzali et all (|2000[) . 

In addition to the decay of 56 Ni as the source of en- 
ergetic 7-photons we include also a number of other ra- 
dioactive decay chains that are important to describe the 
explosion models discussed here. The decays considered 
are: 



Ni - 


■> 57 Co - 


■» 57 Fe, 


Ni - 


+ 56 Co - 


* 56 Fc, 


Fe - 


■> 52 Mn - 


■> 52 Cr, 


Mn - 


+ 51 Cr - 


♦ 51 V, 


Cr - 


■> 49 V - 


■> 49 Ti, 


Cr - 


+ 48 V - 


+ 48 Ti, 


Ti - 


+ 44 Sc - 


+ 44 Ca. 



For the transport of 7-photons we adopt a constant 
gray opacity k 7 = 0.027 cm 2 g _1 and assume that once 
a 7-packagc encounters an interaction it deposits all its 
energy o n the spot to gen erate a package of UVOIR- 
photons (|Swartz et al.l ll995). The amount and rate of 
energy deposit ed in the ejecta is dete r mined foll owing 
the m ethod of lAmbwani fc Sutherland! (|1988l) and iLucvl 
(|2005[) . Some of the decays are accompanied by the emis- 
sion of positrons which are assumed to deposit their ki- 
netic energy in situ. 

The propagation of the resulting uvoiR-photons is 
modeled in a similar Monte Carlo experiment also adopt- 
ing a gray opacity K uvo i r which is, however, parameter- 
ized in terms of the abundance of Fe-group elements in 
the ejecta as 

K uvoiT = 0.25X Fe + 0.025(1 -X Fe ) [cmV 1 ]. (2) 

This paramctrization accounts for the much larger line 
opacity added by the complex ions of Fe-group elements 
compared to lighter elements and has been used in a 
number of previous studies (e.g.. IMazzali et al.l 120011 : 
iMazzali fe Podsiadlowskill200a IShn et al.ll2007t l 

A more detailed description of the modeling procedure 
will be published in an accompanying paper (Saucr et al, 
in prep.) where we study more generally the properties 
of the bolometric light curves from all explosion models 
described in the beginning of this paper. Here we focus 
on the properties of the model which provides the best 
fit to the observed light curve of SN 2005E. 

3.1.1. Results 

We computed bolometric light curves for all explo- 
sion models in Table [TJ Fig. [3] shows the resulting 
light curves for all models in comparison to the quasi- 
bolometric light curve points for SN 2005E (black sym- 
bols) and SN 2002bj (blue symbols). We derived the 
bolometric light curve points for SN 2005E from the 
photometric data published in IPerets et al.l (I2010aj) us- 
ing the procedure developed in lValenti et al.l (|2008l ). In 
addition we made different assumptions to account for 



the unobserved U and JHK bands indicated by the er- 
ror bars on the data points of SN 2005E. For the lower 
limit we took only the observed data from SN 2005E 
while for the upper limi t we added the con tribution of 
the type Ic SN 2007gr (|Hunter et al.ll2009t) . The plot 
symbols correspond to the intermediate light curve ob- 
tained from addin g the U and JHK con tribution of the 
type la SN 2005cf ( |Pastorello et al.ll2007D. Also shown is 
the b olometric light curve of SN 2002b i ([Poznanski et al.l 
120101) . 
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200 400 600 800 
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Fig. 3. — Bolometric light curves of all explosion models in 
comparison to the quasi-bolometric light cu rve points of the ob - 
served supcrnovac SN 2005E (black circles. IPerets et alj|20T0ah . 
and SN 2002bj (blue circles, IPoznanski et all 12010 1 7 . The error 
bars indicate the range of the bolometric light curve points using 
different assumptions for unobserved wavelength bands (see text). 

All models have a comparable rise time to maximum 
of about 7 days. Aside from the absolute brightness at 
maximum light they show strong differences in their late- 
time behavior. The large range of different light curve 
properties of those models is primarily due to the vari- 
ation in the radioactive species synthesized in the ex- 
plosions which have different energy output and decay 
timescales. In addition, the models have slightly different 
density structures resulting from differences in mass and 
kinetic energy. The latter variation affects the timescales 
for the diffusion of trapped photons from the ejecta and 
the brightness at peak. A more in-depth comparison of 
the different light curve models will be presented in a 
forthcoming publication. Here we focus merely on the 
ability of the models to provide a good fit to the ob- 
served data of SN 2005E. The best fit to the data points 
of SN 2005E is provided by the model C0.45HE.2 which 
has a WD mass of 0.45M© and a He-layer of O.2M . The 
model is still somewhat faint at peak and clearly fades 
more rapidly after maximum light than the observed light 
curve of SN 2005E. The rise time of the model light 
curve to maximum light is 6.78 d. Unfortunately the 
pre-maximum data for SN 2005E is sparse and a good 
estimate of the rise time is not available; it seems that 
the rise time in our model is somewhat sho rter than the 
rise tim e of ~11 d estimated for SN 2005E (|Perets et al.l 
12010a!) , but the errorbars on the latter are difficult to es- 
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timate with the available data. None of our models light 
curves show similarities with the bright and ve ry fast 
light curve of SN 2002bj dPoznanski e t alj |2010t or the 
other fast evolving SNe 1885A and 1939B; (jPerets et all 
l2010d ). not shown). The latter bright fast evolving SNe 
arc unlikely to be consistent with the models studied by 
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Fig. 4. — Model light curves from the explosion model C0.45He.2 
(solid black line) in comparison to the quasi-bolometric light curve 
of SN 2005E (shifted to the same rise time as the model). The other 
curves show the contribution of the different radioactive isotopes 
to the total light curve of this model, each labeled with the first 
isotope of the respective /3 + -decay chain. 

Figure [4] shows the light curve of model C0.45HE.2 
alone (solid black line) in comparison to the observed 
data points (shifted to have the same rise time as the 
model). In addition, the contribution of the different 
/3 + decay chains to the total light curves are shown, la- 
beled with the first isotope of each chain. In this partic- 
ular model 44 Ti is actually the most abundant radioac- 
tive species by mass, however, due to its long half life 
of 60 years, this decay only dominates the light curve 
at late times after ^200 d. The dominant contributor 
to the light emission around peak is the decay chain 
48 Cr^ 48 V^ 48 Ti. The time scale of this decay is domi- 
nated by the half life of the second decay of 15.97 d. The 
contribution of the 56 Ni/ 56 Co decays, the most impor- 
tant source of radiative energy in normal, radioactively 
powered supernova light curves, is two orders of magni- 
tude less than the 48 Cr/ 48 V contribution in this model. 

We also tried to check whether the observed light curve 
shape can be matched in a better way with slightly differ- 
ent contributions from the various radioactive isotopes. 
To test this possibility we varied the relative contribu- 
tions of the different decay chains to the total light curve. 
In Fig. [5] we show three different models in addition to 
the original model (solid gray line). In all models the 
density structure is kept the same, only the relative con- 
tributions of the individual light curves shown in Fig. [4] 
are varied. 

Curiously, if one could increase the energy output from 
the 44 Ti decay by a factor of ~50 this model would pro- 
vide an excellent fit also to the late-time light curve 
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Fig. 5. — Tests if variation of the relative contributions of the 
radioactive elements can explain the observed shape of the light 
curve of SN 2005E. In all models only the relative contributions of 
the individual light curves shown in Fig. f4] are varied. The solid 
gray line refers to the original model, model a was obtained by 
increasing the 44 Ti contribution by a factor of 50. For model b 
52 Fe was enhanced by a factor of 30, model c was obtained by only 
using the 56 Ni contribution times a factor of 75 while ignoring all 
other isotopes. 



of SN 2005E (model a in Fig. [5]). However, because 
44 Ti is already the most abundant Fe-group element in 
this model, this increased power cannot be accomplished 
without changing the mass and density structure of the 
entire model, which would inevitably lead to a completely 
different light curve behavior. Nevertheless, this may in- 
dicate that the light emission of SN 2005E is dominated 
by the 44 Ti decay already at times shortly after maxi- 
mum light. A more conclusive assessment would require 
observations out to much later times of several hundred 
days after the explosion where one can safely assume that 
most of the trapped radiation has been released from the 
ejecta and the short-lived isotopes do not contribute sig- 
nificantly to the light emission anymore. 

The pre-maximum data of SN 2005E also indicate that 
this supernova brightened earlier than the model light 
curve. A larger contribution of the 52 Fe light curve by a 
factor of ^30 could explain the early brightening (model 
b in Fig. [5|. However, given the sparse data available for 
SN 2005E a larger ejecta mass leading to a wider light 
c urve cannot b e exclud ed. 

IPerets et al.l (|2010aD estimate an ejected mass of 
0.275M Q and 0.003M Q of 56 Ni from nebular models of 
SN 2005E which, however, do not consider any other 
sources of radiative energy. Model C0.45HE.2 contains 
1.1 x 1CT 4 M Q of 56 Ni. To reproduce the peak of the 
light curve with 56 Ni alone, a factor of ~ 70 more mass 
of this isotope (i.e., ~7.8 x 10~ 3 Af Q ) is needed (model 
c in Fig. [SJ. However, with the given density structure 
56 Ni alone cannot reproduce the decline after maximum. 
The model light curves decay faster than the observation. 

In summary, the decline of the post-maximum light 
curve of SN 2005E indicates that the ejected mass in 
SN 2005E might be somewhat larger than proposed by 
the explosion models. The 44 Ti decay in the model can 
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explain the shape but not the absolute brightness of the 
late time decline. The pre-maximum light curve also 
hints towards a larger ejecta mass, some contribution to 
the early brightening might be explained by the contri- 
bution of the 52 Fe decay chain. 

3.2. Spectral models 

We used our Monte Carlo spectral synthesis code to 
derive a sequence of spectral models from the explo- 
sion models. T h e cod e is bas e d on the d escription by 
IMazzali fc Lucvl (jl993l) ; iLucvl (|1999[ ). and iStehle et all 
( 20051) and has been successfully used in the past for 
the efficie nt interpretation of observed supernova spec - 
tra (e.g., IMazzali et all [20061 [200l ISauer et al.l l2008h . 
The numerical methods are described in more detail in 
the aforementioned publications, here we discuss only the 
most important assumptions as far as they are important 
for the models. 

The code computes a stationary solution for the radia- 
tive transfer through the supernova ejecta using a Monte 
Carlo method with an approximate non-LTE description 
for the atomic level populations. The underlying density 
structure is taken from the hydrodynamic model for ex- 
plosion and is expanded homologously according to the 
time after explosion for each model, assuming that the 
radiation does not alter the hydrodynamic structure of 
the explosion. The radius-dependent composition is ex- 
tracted from the nucleosynthesis results of the explosion 
model. For radioactive species with appreciable abun- 
dance and relevant decay time scales, the conversion to 
the respective daughter elements is taken into account 
according to the epoch. Specifically, the decay chains 
treated explicitly are: 



57 Ni 


->• 


57 Co 


->• 


57 Fe, 


56 Ni 


— >• 


56 Co 


— >• 


56 Fe, 


52p e 


->• 


52 Mn 


->• 


52 Cr 


48 Cr 


— >• 


48y 


— >• 


48 TL 



The solution of the radiative transfer assumes a 
Schustcr-Schwarzschild situation imposing a sharp in- 
ner boundary where all radiation is emitted as a black- 
body. This implies that the 7-photons and positrons 
from radioactive decays below this inner boundary are 
assumed to thermalize, the energy deposition above the 
inner boundary is not taken into account. In supernova 
ejecta it is generally difficult to make a clear choice of 
where to place the inner boundary because the density 
structure of unbound ejecta tends to be very shallow 
such that the location where the ejecta beco me optically 
thick varies strongly with wavelength (e.g., ISauer et al.l 
120061) . For the models discussed here we chose this lo- 
cation iterativcly based on the dilution of the radiation 
field su ch that the dilution fa ctor W becomes close to 
0.5 (see IMazzali fc Lucvlll993l for a discussion of how W 
is determined in our model). With that choice vq corre- 
sponds to the location where the radiation field becomes 
roughly isotropic. 

The approximation of stationarity also implies that 
we cannot determine the luminosity of the model self- 
consistently because the emitted radiative energy at a 
given time has contributions from both the directly de- 
posited energy from radioactive decays and the radiation 
originating from earlier decays which is stored in the op- 
tically thick ejecta and diffuses out as the ejecta expand. 



Therefore, we use the bolomctric luminosity at a given 
epoch from the models of the bolomctric light curves dis- 
cussed in the previous section. 

Once the ejecta start to become diluted enough that 
they are transparent over a wide range of the spectrum, 
the lower-boundary approximation will fail. For the ex- 
plosion models discussed here this generally happens al- 
ready fairly early at about 10 days after maximum be- 
cause the ejected mass of all models is low. In contrast 
to the light curve models which do not depend on the 
photospheric assumption, we cannot compute meaning- 
ful spectra beyond the epoch where that approximation 
breaks down. 

Another complication involves the non-thermal exci- 
tation of He I. He I has high-lying energy levels such 
that the ejecta temperatures are too low to populate 
those levels thermally to give rise to visible absorp- 
tion lines in the spectrum. The He I absorption fea- 
tures seen in SN lb (hyd rogen-deficient su pernovae with 
visible helium lines, see iFilippenkol 119971 for a detailed 
discussion of different types of supernova spectra) are 
caused by non-thermal excitations from fast electrons 
which result from Compton-scatterings of the 7-ph otons 
relea s ed by radioactive decays in the ejecta (e . g. ILucvl 
H991L iKozma fc Franssor][T99l IMazzali fc Lucvlfl998| ) . 

The spectral code we use in this study does not in- 
clude those non-thermal excitations and ionizations self- 
consistently. Therefore, the models will not show signifi- 
cant He I absorption features even though the explosion 
models predict that the He-rich material is well mixed 
into the zones that have a high abundance of radioactive 
isotopes. Nevertheless, assuming that the opacity con- 
tributed by the relatively few He I lines will not substan- 
tially affect the model, the other aspects of the synthetic 
spectra should predict the observable characteristics of 
the explosion models reasonably well. 

3.2.1. Results 

Fig. [75] shows a series of synthetic spectra from the ex- 
plosion model C0.45HE.2 between day 1 and 16 after ex- 
plosion. The flux of each spectrum in this plot has been 
normalized with respect to the maximum flux to make 
the differences in the spectral shape visible. The lower 
two panels show the two spectra of SN 2005E available 
during the photospheric phase. The Jan 16 spectrum cor- 
responds to an epoch of ^3 d before maximum light, the 
Feb 6 spectrum was taken ~18d after maximum. The 
input values for the luminosity and the inner boundary 
velocity vq for all spectral models in this series are shown 
in Fig. [3 

The most prominent ion dominating most of the spec- 
tral shape at all times is Ti II. Also visible are the char- 
acteristic absorption features from Ca II at ^3900 A and 
around 8300 A. The first spectrum in this series is fairly 
red because of the low luminosity on day 1. Day 2 shows 
then a much bluer spectrum which gradually becomes 
redder as time goes on. All spectra show only very lit- 
tle flux blue of the Ca II absorption at ^3900 A. In this 
region the flux is effectively blocked by a dense forest 
of lines from Fe-group elements, mostly Ti II (cf. Fig 
|8|). Comparing the overall shape of the model spectra 
to the two observed spectra one notices that the mod- 
els show most of the observed features, however, the 
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Fig. 6. — Series of synthetic spectra for the explosion model 
C0.45HE.2 (upper panel). The epoch indicated at each spectrum 
is given relative to the explosion date. The absolute flux scale has 
been normalized to the peak of each spectrum to allow for a com- 
parison of the spectral shapes. The (bolometric) maximum of this 
model occurred 6.78 d after explosion. The lower panels show the 
observed spectra of SN 2005E ~3d before and ~18d after maxi- 
mum light. 



strength of absorptions and emission peaks are not well 
reproduced. The spectra before maximum light show 
a prominent peak at ^4000 A which is not as prominent 
in the observed pre- maximum spectrum. In contrast, the 
observed spectrum shows a prominent feature at ^4800 A 
which is partially suppressed by absorption features in 
the model series. 

In the post-maximum phases the model spectra evolve 
faster than the observed spectrum would suggest. This 
is consistent with the light curve model not reproducing 
the observed late-time behavior of SN 2005E. Note that 
the epoch of the Feb 6 spectrum relative to the observed 
maximum of SN 2005E is later than the relative epoch of 
the last model spectrum in the series. The model spec- 
tra after maximum light show strong re-emission features 
between 7500 and 9000 A. Indication for those features 
are present in the observed spectrum, however at very 
different relative strengths. The strong features in the 




4 6 8 10 12 14 16 

epoch [d after explosion] 

Fig. 7. — Input luminosity (upper panel) and velocity of the 
inner boundary (lower panel) for the spectral models as a function 
of time after explosion. 

model spectra originate from strong absorptions and the 
respective re-emission peaks from Ti II lines. 




4000 5000 6000 7000 8000 9000 

rest wavelength [A] 

Fig. 8. — Comparison of the synthetic spectrum from the model 
C0.45He.2 four days after explosion (corresponding to ~3 days be- 
fore maximum light) to the observed spectrum of SN 2005E. The 
flux of the model has been multiplied by a constant to match the 
observed brightness of SN 2005E at this epoch. The labels above 
the spectra indicate the ions that form the main contribution to the 
respective absorption feature. The labels below the spectra indi- 
cate the positions where the stronger He I lines are expected in the 
model spectrum if the model would treat non-thermal excitation 
of He I (see text). 

Fig. [5] shows a detailed comparison of the model spec- 
trum for day 3 and the observed pre-maximum spectrum 
of SN 2005E. The flux of the model spectrum has been 
scaled by an arbitrary factor to match the luminosity 
of the observed spectrum by eye. In this figure the most 
important ions contributing to observed spectral features 
are labeled above the spectra. In particular Ti II lines 
also contribute to most other features in the spectrum. 
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However, for clarity not all are labeled explicitly. Below 
the spectra the positions are indicated where one would 
expect the strongest He I lines to show up if the model 
included non-thermal excitations of this ion. 

The overall shape of the spectrum is reproduced rea- 
sonably well by the model. The line velocities and widths 
arc overall correct indicating that the ejecta velocities 
observed in SN 2005E can be explained by the explosion 
model. The strong peak at ^4700 A seen in the observed 
spectrum is blocked by a number of strong Ti II lines 
in the model. Consequently, we see a strong re-emission 
peak redwards at ~5500A. The model shows Si II lines 
at 4880, 5800 and 6100 A, which arc not or only very 
weakly present in the observed spectrum. Especially the 
absence of the strong Si II 6347, 6371 absorption in the 
observation indicates that the ejecta of SN 2005E con- 
tain only little Si. The total mass of Si in the ejecta of 
model C0.45HE.2 is 1O" 3 M . 

In summary, it appears that the composition of the ex- 
plosion model can reproduce the observed spectrum, per- 
haps the Ti content is somewhat too high in the model, 
however the features from intermediate mass elements 
such as S and Si arc also stronger in the model than in the 
observation. A caveat is that the radioactive isotopes are 
present in the ejecta at almost all velocities. The effects 
of non-thermal ionizations and excitations caused by the 
decays are not accounted for in the radiative transfer 
model. 



4. CONCLUSIONS 

Limited to spherical symmetry, we explored the conse- 
quences of helium detonations in helium layers accreted 
on CO cores for a ra nge of core masses. Recent studies 
by IShen et aLl (|2010l ) explored similar scenarios. Their 
models, however, focused on CO cores of > 0.6 M©, while 
we studied models with CO cores of < 0.6 M©, and ex- 
tended the range of masses to the lowest reasonable limit 
of 0.45 M©. In all cases we a ssumed that carbo n is not ig- 
nited. Our study and that of lShen et al.l (|2010l ) are there- 
fore complementary. We presented here the scenario of 
a CO core of 0.6 M© with 0.2 M© of helium layer. This 
scenario is comparable t o the lowest mass cores scenarios 
explored bv lShen et al.l (|2010l ). Given the different meth- 
ods and simulatio n codes used, our results and those of 
IShen et al.1 (l2010h are in good agreement, showing quite 
similar light curves and nucleosynthesis products. 

Our results show that below a CO core mass of 0.6 M©, 
the abundance of 56 Ni, produced by the detonation, 
drops rapidly, while those of 40 Ca, 44 Ti, and 48 Cr grow 
rapidly to a few percent. Moreover, the fraction of un- 
burnt helium also grows with decreasing CO core mass. 
These findings are encouraging for the discussed scenario 
to be a viable model for SN 2005E-like objects, and sug- 
gest that the latter, low mass CO core scenarios on which 
we focused are the more likely progenitors of these ex- 
plosions. It is important to note that, in a narrow mass 
range of models, a highly diverse outcome is expected, 
both in terms of light curves and nucleosynthesis. Hence, 



one cannot expect an accurate fit of a model to a specific 
observed supernova. Nevertheless, more detailed calcu- 
lations of the spectra, including the helium lines and late 
nebular spectra are yet to be done. 

Besides the spectrum, the main remaining question 
concerns the multi dimensional aspect of the problem, 
which has a strong impact on whether a consecutive 
CO detonation may follow the helium detonation or not. 
We speculate that in the 2D case and with low CO 
core masses (0.5 M© and under), the density of the core 
around the possible ignition sites will be too low for CO 
ignition. The late evolution to thermonuclear runaway 
and the onset of helium detonation are also issues which 
require at least two dimensional tools. This is due to 
the convective nature of the reactive flow at the base of 
the helium layer, prior to runaway. Dredge up of heavier 
nuclei into the helium layer could in principle alter the 
conditions for a transition from quasi-static burning to a 
detonation. 

In terms of nucleosynthetic products, the large abun- 
dances of 40 Ca and 44 Ti produced in our detailed ex- 
plosion models confirm similar results obtained with the 
si mplified one-z o ne nucl eosynthetic calculation discussed 
in iPerets et al.l ()2010al ). As discussed there, our de- 
tailed model results could therefore have important im- 
plications for the calcium enrichment of the interstellar 
medium as well as the the production of positrons and 
the production of the 511 keV annihilation line. 

Finally, if SN 2005E type objects are produced by he- 
lium detonations on small CO cores, there should also 
be a subclass of objects formed by helium detonations 
on larger cores. In that class the spectrum should be 
dominated by lines of radioactive nickel, and the proba- 
ble ignition of the CO core should produce much more 
energetic and metal rich objects. We therefore encour- 
age observers to search for such supernovae in the near 
future. 
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TABLE 1 

Parameters of the Simulated Initial Configurations 



Model 


M co 




P6c 


T 7c 


P6He 




C0.45HE.2 


0.45 


0.2 


3.81 


1. 


0.543 


20. 


CO.5HE.15 


0.5 


0.15 


3.92 


1. 


0.402 


20. 


C0.5HE.2 


0.5 


0.2 


5.06 


1. 


0.678 


20. 


C0.5HE.3 


0.5 


0.3 


8.50 


1. 


1.391 


20. 


C0.55HE.2 


0.55 


0.2 


6.72 


1. 


0.845 


20. 


CQ.6HE.2 


0.6 


0.2 


8.81 


1. 


1.032 


20. 



Note. — mass units— Mq, pe c is the central density in units of 10 Q g cm 3 , pe,He is the density at the base of the helium 
shell in same units, Tj c and Tjhe are the corresponding temperatures in units of 10 7 K . 
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TABLE 2 

Explosion Energy and Nucleosynthesis Products 



Isotope 


C0.45HE.2 


CO.5HE.15 


C0.5HE.2 


C0.5HE.3 


C0.55HE.2 


C0.6HE.2 


Ek 


0.178 


0.096 


0.201 


0.460 


0.226 


0.242 


4 He 


1.1E-01 


9.8E-02 


1.0E-01 


9.7E-02 


9.3E-02 


8 3F-02 

O • U J— J Uii 


20 Ne 


q 8E-0fi 

Zj . O 1 ' VJ U 


1.2E-05 


6 ^iF-Ofi 

V) . J J_J U U 


3 9E-0R 

■ ) . ij J_J U U 


4.7E-06 


3.7E-06 


23 Na 


1.2E-09 


2.0E-09 


6.0E-10 


1.4E-09 


3.0E-10 


1.6E-10 


24 Mfr 


4.4E-05 


fi 2F-0 C > 

U ■ £i 1 J U J 


2.7E-05 


1.3E-05 


1.8E-05 


1.4E-05 


27 Al 


9.9E-07 


1.1E-06 


7.0E-07 


3.1E-07 


4.5E-07 


4.0E-07 


28 Si 


1.0E-03 


1.5E-03 


^ 8F 04 


1.7E-04 


1 SF 04 


2.7E-04 


31p 


1.8E-05 


2 2F-0 c i 


1.2E-05 


2 fiE-O 1 ; 

__ . 1 1 1 J U J 


8 3E-06 

O . JiJ UU 


fi fiF-Ofi 

U . UJJ UU 


32 g 


4.6E-03 


fi 9E-03 


2.7E-03 


5.5E-04 


1.7E-03 


1.1E-03 


35 CI 


1.4E-04 


1.3E-04 


1.5E-04 


1.3E-04 


1.2E-04 


fi fiF-O^i 

U . UJJ UJ 


36 Ar 


J . J i ' u i-i 


9.1E-03 


3 2E-03 

J . 1 J U J 


6.0E-04 


1.9E-03 


1.4E-03 


39 K 


Q 4F OA 


fi &F ru 

U . OJ_J Url: 


1.1E-03 


4 QF 04 


8.2E-04 


4 1 F-fU 

■rt. J_ J_J U1 


40 Ca 


3.4E-02 


2.9E-02 


2.1E-02 


^ c iF-03 

J . J 1 ' VJ J 


1.3E-02 


9 OF-03 

•J . UJJ U J 


41 Ca 


3 3E-0fi 


1.8E-06 


5.4E-06 


1.3E-06 


4.9E-06 


2.3E-06 


42 Ca 


J . J J_J U U 


LIE- 06 


1.8E-05 


8 3F-0R 


1.5E-05 


4.7E-06 


43 Ca 


8.4E-05 


2 3F-0 c i 

_■ i J J_ J U J 


1.1E-04 


J . J J_J U J 


8 OF-O^i 

(J . U1J V) ' ) 


^ 3F-0 c i 

J . JJ-J U J 


44 Ca 


1.4E-06 


1.7E-07 


1.3E-06 


2.3E-07 


7.9E-07 


5.2E-07 


45 Sc 


1.6E-05 


1 8F Ofi 

J. .OJ_J uu 


1.5E-05 


2.1E-06 


K 9F Ofi 

O. UU 


1 SF Ofi 




3 3E-02 


3 2E-03 

J . £i J_ J UO 


3.1E-02 


^ 8F-03 

J . O J_J U J 


2 OF-02 

_ . U J—J Ui 


1.3E-02 




4.0E-06 


5.8E-07 


4.2E-06 


7.8E-04 


2 8F-06 

_ . O J—J UU 


fi QF-O^i 

U . C/J-J U J 


47 Ti 


3.7E-04 


1.8E-05 


7.4E-04 


1.8E-03 


6.0E-04 


4.6E-04 


48 Ti 


9.4E-05 


8.0E-07 


3.5E-04 


8 8E-0 1 ; 

O . O J_J U J 


2.7E-04 


2.2E-04 


49<pj 


1.4E-07 


1.5E-09 


4 fiF 07 


1.3E-07 


2.6E-07 


1.2E-07 


47y 


1.1E-06 


9 OF n& 


^ 2F Ofi 


fi i F ns 

U. J_J_j UO 


Q 9F flfi 

7J . ZJ-j UU 


9 QF Ofi 

£l . ,7 J_J UU 


48y 


3 2E-(13 


2.7E-05 


1.2E-02 


3 0E-03 


9.4E-03 


7.4E-03 


49 V 


U . UU U J 


6.2E-07 


2.0E-04 


J. J1j UJ 


1.1E-04 


J i __ J— J U J 


51y 


1.0E-05 


Fi QE-08 


4.7E-05 


fi ^F-O'i 

V) . U1j UJ 


3 QE-O^i 


2.7E-05 


48 Cr 


2 3F-03 


2 OF-O^ 

_■ * U J— J U J 


8.7E-03 


2 2F-03 

_ , ilj UJ 


7.0E-03 


J . JJ-J U J 


49 Cr 


1.4E-07 


4 4F 1 

'-1 . 4l1_J J_U 


7.7E-07 


3.1E-08 


9 9F Dfi 

_ . UU 


1 4F Dfi 

J_ . tJjj uu 


50 Cr 


2 flFj-Ofi 

_ . Ul-J uu 


2.3E-07 


6.1E-05 


3.8E-04 


3 6E-0 t S 

J.UJJ UJ 


3 2F-0 c i 

J > J_ J U J 


51 Cr 


3.6E-04 


2. IE- 06 


1 7E-03 


2 3E-03 


1.4E-03 


9.4E-04 


51 Mn 


2.4E-07 


2.2E-10 


2 ^iF-Ofi 

— . J 1 J U U 


2.3E-07 


1.6E-05 


1.4E-05 


52 Mn 


7.8E-04 


2 0E-06 

_ • U J— J U I. ' 


9.1E-03 


1.6E-02 


2.4E-02 


2.1E-02 


^ 3 Mn 


c OF 0^ 


1 2E-07 


K Oil 
J. ( JZj- U '-i 


9 fiP OA 


1 2E-03 


4 7F 04 


51 Fe 


9.7E-08 


4.0E-11 


7.9E-07 


2.2E-07 


2.5E-06 


2.3E-06 


52 Fe 


8.5E-05 


2.1E-07 


9.9E-04 


1.7E-03 


2.7E-03 


2.6E-03 


53 Fe 


2.6E-08 


2.9E-12 


3.3E-07 


7.7E-09 


5.7E-06 


1.6E-05 


54 Fe 


1.7E-05 


4.2E-08 


1.5E-04 


3.7E-04 


2.7E-04 


1.2E-04 


55 Fe 


7.8E-05 


1.2E-07 


7.7E-04 


1.8E-03 


2.4E-03 


1.5E-03 


56 Fe 


9.5E-08 


1.0E-10 


1.0E-06 


7.9E-05 


1.1E-05 


2.7E-05 


57 Fe 


2.2E-08 


2.8E-11 


1.7E-07 


1.0E-05 


1.4E-06 


3.1E-06 


55 Co 


3.9E-05 


5.8E-08 


3.8E-04 


9.1E-04 


1.2E-03 


7.8E-04 


56 Co 


1.5E-05 


1.5E-08 


1.8E-04 


1.5E-02 


2.0E-03 


5.2E-03 


57 Co 


1.3E-05 


1.7E-08 


1.1E-04 


6.4E-03 


8.3E-04 


2.0E-03 


56 Ni 


9.9E-05 


9.8E-08 


1.2E-03 


1.1E-01 


1.4E-02 


3.7E-02 


57 N ; 


1.8E-05 


2.3E-08 


1.5E-04 


9.1E-03 


1.2E-03 


2.8E-03 


58 Ni 


1.6E-05 


3.0E-08 


1.4E-04 


1.8E-02 


8.7E-04 


2.1E-03 


59 N ; 


1.4E-07 


4.4E-10 


1.8E-06 


3.4E-05 


8.6E-06 


1.8E-05 


60 N ; 


4.0E-09 


2.1E-11 


7.7E-08 


6.0E-07 


4.6E-07 


1.0E-06 



NOTE. — The products listed include the He layer only, and are given at the beginning of radiative transfer calculation 
(t — 10 5 s). Between He and K, only the most abundant isotope is included, from Ca and up all isotopes exceeding 1O -6 M0 
in at least one model are listed. EK - final kinetic energy in units of 10 51 ergs, isotopes - in solar mass. 
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TABLE 3 



End product as a function of initial 12 C mass fraction 


Xq (initial) 


N 


End-product ( 12 C + No) 


0.6 


2 


20 Ne 


0.5 


3 


24 Mg 


0.429 


4 


28 Si 


0.375 


5 


32 g 


0.33 


6 


36 Ar 


0.30 


7 


40 Ca 


0.273 


8 


44 Ti 


0.25 


9 


48 Cr 


0.231 


10 


52 Fe 


0.214 


11 


56 Ni 



